Equation of State of Gluon Plasma from Fundamental Modular Region 
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Despite considerable practical success in dealing with the gluon plasma, finite-temperature perturbation the- 
ory suffers at the fundamental level from infrared divergences discovered by Linde. However if gauge or Gribov 
copies are properly eliminated from the physical state space, infrared modes are strongly suppressed. We de- 
scribe the gluon plasma in zeroth order as a gas of free quasi-particles with a temperature-independent dispersion 

relation of Gribov type, E{k) = ^Jk^ + that results from the reduction of the physical state space. The ef- 
fective mass ^ controls infrared divergences and allows finite calculable corrections. The equation of state of 
this gas is calculated and compared with numerical lattice data. 



I. INTRODUCTION 



The equation of state (EOS) of the quark-gluon plasma is 
of vital interest at present because of the exciting RHIC ex- 
periments at Brookhaven National Laboratory [1]. Here we 
shall be concerned with the pure gluon plasma, leaving quarks 
for another occasion. The EOS of the gluon plasma is rather 
precisely known from numerical studies in lattice gauge the- 
ory [2], and there exist excellent phenomenological fits to 
the EOS of the gluon plasma in the deconfined phase [3]. 
In contrast with this practical success, it was discovered by 
Linde [4] as far back as 1980 that finite-temperature pertur- 
bative QCD suffers at a fundamental level from infrared di- 
vergences, which suggests that finite-temperature perturba- 
tion theory neglects an essential feature of QCD. It has been 
proposed [5] to control these divergences by introducing a 
magnetic mass m ^ g^T. It has also been proposed [6] that 
the divergences stem from an inadequate application of the 
principle of gauge equivalence. Indeed in 1978, two years 
before Linde's discovery of the infrared divergence [4], Gri- 
bov showed that infrared modes are strongly suppressed when 
gauge equivalence is imposed at the non-perturbative level [7] . 

According to the principle of gauge equivalence, two differ- 
ent configurations that are related by a gauge transformation, 
A2 — ^Ai, represent the same physical state. The gauge trans- 
formation, ^A, = g^^Aig + g^^djg, is highly non-linear, and 
the physical state space (which is the quotient of the space 
of connections modulo gauge transformations) is not a linear 
space. This geometric property of gauge theory shows up in 
practice when one fixes a gauge, which is simply a coordi- 
natization of the physical state space. Two or more different 
configurations may be gauge equivalent even though both sat- 
isfy the linear Coulomb gauge condition, ^^Lj 3,A/ = 0. When 
enumerating physical states, only one of these "Gribov" or 
gauge copies should be counted, so the space of physical states 
is reduced to the fundamental modular region (FMR), a region 
that is free of Gribov copies. Gribov [7] found that the disper- 
sion relation E{k) — k gets modified because of the reduction 



of the physical state space, and he obtained instead 



(1.1) 
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where k = |k|, and M is a QCD mass scale. The reduction 
of the physical state space was originally proposed as an es- 
sential feature of the confinement mechanism [7, 8]. However 
statistical mechanics is primarily a matter of counting states, 
and the reduction of the physical state space required by the 
gauge principle influences the EOS at all temperatures. Here 
we shall be concerned with its effect in the deconfined phase. 

It is known from numerical studies [2] that at high tem- 
perature the EOS of the gluon plasma approaches the Stefan- 
Boltzmann law, e = 3p = 3csbT'^, where e is the energy per 
unit volume, p is the pressure, T is the temperature, and 
csB = ^{N^ — 1) in SU(N) gauge theory. Thus it seems rea- 
sonable to describe the gluon plasma at high temperature in 
first approximation as a gas of non-interacting quasi-particles. 
We shall describe the quasi-particles by the Gribov disper- 
sion relation (1.1) or a similar one, for E{k) is only approx- 
imately known. Fortunately the results obtained hold under 
rather general conditions on E{k). We call the gas of non- 
interacting quasi-particles, with modified dispersion relation, 
the FMR gas. 

The FMR gas closely resembles early phenomenologi- 
cal models of the gluon plasma [9], characterized by a 
temperature-independent infrared cut-off K, because the ef- 
fective mass which is large at low k, provides an infrared 
cut-off at k ^ M. Those models, which were developed ap- 
proximately a decade after Gribov's paper, were inspired by 
the idea that low-momentum gluons are confined but high- 
momentum gluons are effectively free, and they successfully 
capture the gross features of the EOS of the gluon plasma. 
Phenomenological models were subsequently improved [3], 
by making the cut-off K{T) temperature dependent and in 
other ways, to give a good fit in the transition region, at the 
cost of additional parameters. The FMR gas is not intended to 
compete in precision with these improved phenomenological 
models. Rather it is conceived as the zeroth order approxi- 
mation that allows calculable corrections by an iterative pro- 
cedure, for example of Dyson-Schwinger type (presently be- 
ing developed). Successive iterations yield perturbative-type 
corrections to the FMR gas that are expected to give good 
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results at high T . They are similar to thermal perturbation 
theory [10, 11], but the infrared divergences of thermal per- 
turbation theory [4], that may be controlled by introducing a 
magnetic mass m ^ g^T [5], are here controlled by the effec- 
tive mass It has been noted by Rischke et al [9] that in- 
frared divergences are also controlled when the cut-off model 
is iterated, and that non-analytic corrections of order are 
absent. 

Many aspects of the Gribov scenario have been verified in 
recent investigations. Infrared suppression of the gluon propa- 
gator in Coulomb gauge has been observed in numerical sim- 
ulation [12], but less strongly in [13], and in Landau gauge 
in 3-dimensions [14]. It has also been found in variational 
calculations in Coulomb gauge [15], and in Schwinger-Dyson 
calculations in Coulomb [16| and in Landau [17] gauge. A 
long-range color-Coulomb potential was found in numerical 
simulations in Coulomb gauge in the deconfined phase [ 1 8 1 . 

The questions we wish to address here are: (i) What is the 
equation of state of the FMR gas? (ii) How does it compare 
with the EOS that is known from numerical studies? 



integration by parts, one has 

^ 37t2 exp[p£(A:)] - 1 dk ' ^ ' 

From the last two equations, we obtain for the trace anomaly 

= e-3;?=(G2)o-(G2)r, 

^~^^Jo "'^exp[p£(/t)]-l dk\~^)' ^^-^^ 

where = (2^^)^^(— p)G",vG"^" is the gluon conden- 
sate [19]. The integrand is positive when ^ (^^^^ < 0> which 
is a sufficient condition for e{T) > 3p{T). We also have 
e = r^^(^), from which it follows that the trace anomaly 
may be written Q = e—3p = T^-^{^). Upon integration this 
yields 

poo 

P = csrT^-T'^ dT' T'-^ e{T'). (2.6) 



II. EOS OF THE FMR GAS 

In the Stefan-Boltzmann limit, the degrees of freedom for 
each gluon momentum k are the two states of polarization, 
each with color multiplicity (A'^ — 1)- These are precisely the 
degrees of freedom in Coulomb gauge, and we shall use this 
gauge for our calculation. It is a "physical" gauge without 
negative metric states, and all constraints are satisfied identi- 
cally. Although the Coulomb gauge is not manifestly Lorentz 
covariant, this is not necessarily a disadvantage in the de- 
confined phase because at finite T the heat bath provides a 
preferred Lorentz frame, and the manifest symmetries of the 
Coulomb gauge are the synmietries of the physical problem at 
hand. 

The partition function of the FMR gas is given by the 
Planck distribution, 

Z = n[l-exp(-P£„)]-', (2.1) 

n 

where (3 = l/T is the inverse temperature, and n = (k, A-,a), 
where k is 3-momentum, X = 1,2 is polarization, and a = 
1,...{N^ - 1) is color. With l,k,X,a ^ (N^ - l)V%-^ J dkk^, 
where V is the 3-volume, this gives 

InZ = -{N^ - iWn-^ ( dkk^ ln{ 1 - exp[-^E{k)] }, 
Jo 

(2.2) 

and, with energy density e = — 7 we have 

(N^-\) r k'^ 
e = ^1—^ / dk — , E{k). (2.3) 

For a homogeneous system, the pressure is given hy p = 
= y InZ = — /, and the entropy density hy Ts ~ e + p, 
where f = F /V is the free energy per unit volume. After an 



III. FMR GAS AT HIGH TEMPERATURE 

Suppose that the leading deviation of E{k) from k at high k 
is expressed by a power law, 

E{k)/k=l+c/0+... . (3.1) 

For the Gribov dispersion relation E{k) = i/P'+M^T^ 
has y = 4, whereas for a gluon mass, E{k) = vP+ff?, one 
has y = 2. The gluon condensate has dimension 4, which leads 
one to expect y = 4, whereas if there were a condensate of di- 
mension 2, one would expect J =2. The asymptotic behavior 
of the EOS is qualitatively different for y greater or less than 3. 
We suppose y > 3, and consequently the deviation from the ul- 
traviolet behavior £■ (A:) = kis soft. 

For y > 3, the asymptotic high-r limit of the trace anomaly 
is obtained from the substitution exp[^p^] — 1 — > in (2.5). 
This substitution cannot be made in the integrals for p and e 
because they would diverge. It gives a Hnear asymptotic trace 
anomaly, 

e=L7' -1-0(1), (3.2) 

where 

7t2 Jo E{k) dk\ k J 

poo 

= 3(N^-l)%-^ dkk^ ln[E(k)/k], (3.3) 
Jo 

is an integral that converges for y > 3 (as we have supposed), 
and is positive for E{k) > k. These are sufficient conditions 
for the hnear asymptotic form (3.2) with L > 0. Although L 
describes the high-T limit of the FMR gas, the last integral 
involves all momenta k. For the special case of the Gribov 
dispersion relation (1.1), one obtains 

L={N^-l){nV2)-^ m\ (3.4) 
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which is proportional to although only appears in the 
dispersion relation. 

The pressure at high T is obtained from (2.6) which yields 

P = csbT'^-{1/3)LT + 0{1), (3.5) 

Thus for the FMR gas, the leading deviation of the pressure 
from the Stefan-Boltzmann law is linear in T. However this 
linear term — and only a linear term — is annihilated in the 
formula for the energy density 8 = r^^(^), which gives 

e = 3csBT^ + 0{l). (3.6) 

An EOS of this type was obtained as a fit to the lattice data at 
high temperature in [20] . For the speed of sound one obtains 

IV. COMPARISON WITH NUMERICAL EOS 

The EOS of the FMR gas at high T is not sensitive to the 
exact form of E{k) because (3.5) and (3.6) hold as long as 
Y > 3 and E{k) > k. To compare with the numerical data, we 
take the Gribov dispersion relation (1.1). The unknown mass 
scale M is determined by fitting the anomaly = e — 3/:' at 
high T, because the corrections to this quantity are expected 
to be small. [Indeed, taking thermal perturbation theory as 
a guide [10], we note that the leading correction to p is of 
order Ap ~ g^{T)T'^ ^ The anomaly is given by 9 = 
r^^(^), so the corresponding correction to the anomaly, 
AO = J.(f£), is of order T'^{^) = -^- g\T)T\] 
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FIG. 1: Numerical and analytic plot of (e — 3p) /T'^. 
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FIG. 2: Numerical and analytic plots ofe/T*, 3/4 s/T^, 3p/T'^. 



The numerical data of [2] for ^ are represented by the 
black interpolating curves in Fig. 1. (More recent studies in- 
clude dynamical fermions or a chemical potential that cannot 
be described by the FMR gas.) The data for Ni — 6 and 8 
agree, and were interpreted as continuum values [2]. The red 
dots are obtained from the analytic formula (2.5), with mass 
scale set at M = 2.6rc- by fitting at high T, where Tc is the 
transition temperature. The relatively large difference in the 
transition region between the FMR gas and the numerical data 
for Q = e—3p occurs because the FMR gas does not exhibit 
a sharp phase transition, whereas for pure SU(3) gauge theory 
there is a first order phase change, so e is discontinuous while 
p is continuous. We do not attempt to estimate the error of 
M because perturbative-type corrections to the FMR gas have 
been neglected. From the value — 0.625 \/o of [2], where a 
is the string tension, one gets M — 1 or M = 705 MeV, 
where the string tension for the quarkless theory is defined to 
be = 440 MeV. 

The numerical data of [2] for ^ and |§ are are dis- 
played as black interpolating curves in Fig. 2, where jj = 

^j^. The horizontal line represents the Stefan-Boltzmann 
EOS. For the FMR gas, |§ and ^ approach the Stefan- 
Boltzmann limit like jj and ^ respectively, being quite 
close to it at the highest temperature displayed, T = 5Tc, 
whereas the gluon plasma approaches the Stefan-Boltzmann 
limit more slowly. The difference between the FMR gas and 
the gluon plasma in the range 2Tc to 571 appears attributable 
to perturbative-type corrections of moderate size. In standard 
thermal perturbation theory [10] these are of leading order 
g'^{T) ~ jj^, but diverge at order g^, whereas corrections to 
the FMR gas are expected to be calculable. 



V. DISCUSSION 

Above the transition region, the EOS of the FMR gas gives 
a good description of the most prominent feature of the gluon 
plasma which is the rapid drop of the pressure compared to the 
energy from the Stefan-Boltzmann value, as T decreases from 
infinity. The linear asymptotic trace anomaly (3.2) provides 
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a ready explanation for this, that holds for any quasi-particle 
model with y > 3 and ^ (^^y^^ < although other fits are 
certainly not excluded. The FMR gas is not exact even at high 
temperature because of perturbative-type corrections. We ex- 
pect however that they are calculable and of moderate size 
above the transition region. 

The transition region is not so well described by the FMR 
gas. There is no sharp phase transition because p{T) is an 
analytic function. Moreover the dependence on N [of SU(N)] 
is only through the coefficient (A^^ — 1)> whereas even the or- 
der of the phase change depends on A', being second order for 
SU(2) and first order for SU(3). An analysis of the phase tran- 
sition based on center symmetry is given in [21]. However 

the dispersion relation E{k) = \Jk^ + ^ has a minimum en- 
ergy iiniin = V^M, SO the thermodynamic functions e, p and s 
of the FMR gas are exponentially small for T < VlM = 997 
MeV (forM = 705 MeV). The mass of the lightest glueball is 
of order 1 GeV, so the thermodynamic functions of the FMR 
gas are exponentially small where they are supposed to be. 
Altogether the FMR gas, with a single parameter M which is 
the mass scale, is competitive with cut-off phenomenological 
models [9]. 

It may be helpful to compare the FMR gas with the more 
perfected phenomenological models [3]. The dispersion re- 
lation, E{k), which reflects the restriction to the FMR, is T- 
independent, and consequentiy the thermodynamic identities 
are automaticaUy satisfied using Planck's formula. By con- 
trast the phenomenological dispersion relations E{k, T) of [3] 



are T-dependent. When used in Planck's formula, the T- 
dependent dispersion relation requires the introduction of an 
additional "background field" B{T) to assure that the thermo- 
dynamic identities are satisfied. The improved phenomeno- 
logical models of [3] are more precise than the FMR gas es- 
pecially in the transition region, but they also require several 
parameters. We have fit only the QCD mass scale M. 

However the FMR gas is not intended to be a precise phe- 
nomenological model, but rather to provide a useful starting 
point, well founded in the principles of gauge theory, that al- 
lows calculable, moderate size corrections at high T . Athough 
it is defined by Gribov's dispersion relation of 1978, the FMR 
gas has two important properties that were later independently 
reinvented: (i) Its EOS closely resembles simple phenomeno- 
logical models of the gluon plasma [9]. (ii) The effective 
mass ^ controls infrared divergences in higher order correc- 
tions so it is not necessary to introduce the magnetic mass 
M'^ g-^T [5] for this purpose. 
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